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Abstract 



The compelling experimental evidences for oscillations of solar, reactor, atmospheric and accelerator 
neutrinos imply the existence of 3- neutrino mixing in the weak charged lepton current. The current 
data on the 3-neutrino mixing parameters are summarised and the phenomenology of 3-u mixing 
is reviewed. The properties of massive Majorana neutrinos and of their various possible couplings 
are discussed in detail. Two models of neutrino mass generation with massive Majorana neutrinos 
- the type I see-saw and the Higgs triplet model, are briefly reviewed. The problem of determining 
the nature - Dirac or Majorana, of massive neutrinos is considered. The predictions for the effective 
Majorana mass | ( m ) | in neutrinoless double beta ((/3/3)o^ — ) decay in the case of 3-neutrino mixing and 
massive Majorana neutrinos are summarised. The physics potential of the experiments, searching for 
(PP)ov— decay for providing information on the type of the neutrino mass spectrum, on the absolute 
scale of neutrino masses and on the Majorana CP-violation phases in the PMNS neutrino mixing 
matrix, is also briefly discussed. The opened questions and the main goals of future research in the 
field of neutrino physics are outlined. 



1 Invited review article for the special issue on Neutrino Physics of the journal Advances in High Energy Physics. 
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1.1 Introduction: The Three Neutrino Mixing (an Overview) 

It is a well-established experimental fact that the neutrinos and antineutrinos which take part in the 
standard charged current (CC) and neutral current (NC) weak interaction are of three varieties (types) 
or flavours: electron, v e and v e , muon, i>„ and z7„, and tauon, v T and v T . The notion of neutrino type 
or flavour is dynamical: v e is the neutrino which is produced with e + , or produces an e~ in CC weak 
interaction processes; Uu is the neutrino which is produced with or produces fi~, etc. The flavour 
of a given neutrino is Lorentz invariant. Among the three different flavour neutrinos and antineutrinos, 
no two are identical. Correspondingly, the states which describe different flavour neutrinos must be 
orthogonal (within the precision of the current data): {vvpi) = 5m, (Pi'\vi) = Sin, = 0. 

It is also well-known from the existing data (all neutrino experiments were done so far with relativis- 
tic neutrinos or antineutrinos), that the flavour neutrinos v\ (antineutrinos vi), are always produced 
in weak interaction processes in a state that is predominantly left-handed (LH) (right-handed (RH)). 
To account for this fact, v\ and v\ are described in the Standard Model (SM) by a chiral LH flavour 
neutrino field vn(x), I = e,fi,T. For massless vi, the state of v% (pi), which the field vil(x) annihilates 
(creates), is with helicity (-1/2) (helicity +1/2). If v\ has a non-zero mass m(yi), the state of vi (pi) 
is a linear superposition of the helicity (-1/2) and (+1/2) states, but the helicity +1/2 state (helicity 
(-1/2) state) enters into the superposition with a coefficient oc m{yi)jE, E being the neutrino energy, 
and thus is strongly suppressed. Together with the LH charged lepton field li(x), vil(x) forms an 
SU{2)l doublet in the Standard Model. In the absence of neutrino mixing and zero neutrino masses, 
vil(x) and Il(x) can be assigned one unit of the additive lepton charge L[ and the three charges L\, 
I = e, //, t, are conserved by the weak interaction. 

At present there is no compelling evidence for the existence of states of relativistic neutrinos 
(antineutrinos), which are predominantly right-handed, vr (left-handed, vi). If RH neutrinos and LH 
antineutrinos exist, their interaction with matter should be much weaker than the weak interaction 
of the flavour LH neutrinos v\ and RH antineutrinos v\, i.e., vr (pi) should be "sterile" or "inert" 
neutrinos (antineutrinos) \X\. In the formalism of the Standard Model, the sterile vr and vi can be 
described by SU(2)i singlet RH neutrino fields vr(x). In this case, vr and i>i will have no gauge 
interactions, i.e., will not couple to the weak and Z° bosons. If present in an extension of the 
Standard Model (even in the minimal one) , the RH neutrinos can play a crucial role i) in the generation 
of neutrino masses and mixing, ii) in understanding the remarkable disparity between the magnitudes 
of neutrino masses and the masses of the charged leptons and quarks, and iii) in the generation of 
the observed matter-antimatter asymmetry of the Universe (via the leptogenesis mechanism [2], see 
also, e.g., [3]). In this scenario which is based on the see-saw theory [3], there is a link between the 
generation of neutrino masses and the generation of the baryon asymmetry of the Universe. The 
simplest hypothesis (based on symmetry considerations) is that to each LH flavour neutrino field 
vil(x) there corresponds a RH neutrino field vir(x), I = e,fi,r, although schemes with less (more) 
than three RH neutrinos are also being considered (see, e.g., [5]). 

The experiments with solar, atmospheric, reactor and accelerator neutrinos (see [6] and the refer- 
ences quoted therein) have provided compelling evidences for flavour neutrino oscillations [H IE] - 
transitions in flight between the different flavour neutrinos v e , v^, v T (antineutrinos v e , v^, v T ), caused 
by nonzero neutrino masses and neutrino mixing. As a consequence of the results of these experiments 
the existence of oscillations of the solar v e , atmospheric and v^, accelerator (at L ~ 250 km, 
L ~ 730 km and L ~ 295 km, L being the distance traveled by the neutrinos) and reactor v e (at 
L ~ 180 km and L ~ 1 km), was firmly established. The data imply the presence of neutrino mixing 
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in the weak charged lepton current: 



n 



Ccc = ~ ~ E l ^ x ) Wil(x) W a \x) + h.c. , u lL (x) = £ ^i^L(x), (1.1) 



where vil{x) are the flavour neutrino fields, z^l(30 is the left-handed (LH) component of the field of 
the neutrino Uj having a mass mj, and U is a unitary matrix - the Pontecorvo-Maki-Nakagawa-Sakata 
(PMNS) neutrino mixing matrix [HOE], U = Upmns- All compelling neutrino oscillation data can 
be described assuming 3-neutrino mixing in vacuum, n = 3. The number of massive neutrinos n can, 
in general, be bigger than 3 if, e.g. there exist right-handed (RH) sterile neutrinos and they mix with 
the LH flavour neutrinos. It follows from the current data that at least 3 of the neutrinos Vj, say v\, 
V2-, ^3, must be light, m.1,2,3 ^ 1 eV, and must have different masses, mi ^ m<i ^ m^. At present there 
are no compelling experimental evidence for the existence of more than 3 light neutrinos. Certain 
neutrino oscillation data exhibit anomalies that could be interpreted as being due to the existence 
of one or two additional (sterile) neutrinos with mass in the eV range, which have a relatively small 
mixing ~ 0.1 with the active flavour neutrinos (see, e.g., [9] and the references quoted therein.) 

In the case of 3 light neutrinos we will concentrate on in this Section, the neutrino mixing matrix 
U can be parametrised by 3 angles and, depending on whether the massive neutrinos Uj are Dirac or 
Majorana PH] particles, by 1 or 3 CP violation (CPV) phases [TT1 [T2l [TB] : 



^21 .■ ".■{! 



U = VP, P = diag(l,e*— ,e l ~ ), (1.2) 

where 021 and 031 are the two Majorana CPV phases and V is a CKM-like matrix containing the 
Dirac CPV phase 8, 



V 



I C12C13 S12C13 si2,e~ %& \ 

-S12C23 - Cl2S 2 3Sl3e* 5 C12C23 - Sl2S23Sl3e*' 5 S23C13 
V S12S23 - Ci 2 C 2 3.Sl3e i<5 -C12S23 - Sl2C23Sl3e i<5 C23C13 / 



1.3) 



In eq. (1.3), c%j = cosOij, Sij = sinOij, the angles Oij = [0,7r/2], 5 = [0,2ir] and, in general, < 
< 2tt, j = 2,3 [H]. If CP invariance holds, we have 5 = 0,ir, and [13 [13], a 2 i(3i) = k^ir, 
feCJ = 0,1,2,3,4. 

Thus, in the case of massive Dirac neutrinos, the neutrino mixing matrix U is similar, in what 
concerns the number of mixing angles and CPV phases, to the CKM quark mixing matrix. The 
presence of two additional physical CPV phases in U if Vj are Majorana particles is a consequence of 
the special properties of the latter (see, e.g. refs. [HI IS]). On the basis of the existing neutrino data 
it is impossible to determine whether the massive neutrinos are Dirac or Majorana fermions. 

The neutrino oscillation probabilities depend, in general, on the neutrino energy, E, the source- 
detector distance L, on the elements of U and, for relativistic neutrinos used in all neutrino experiments 
performed so far, on the neutrino mass squared differences Am| = (mf — m|), i ^ j (see, e.g., 
ref. [H [IE]). In the case of 3-neutrino mixing there are only two independent neutrino mass squared 
differences, say Am|i ^ and Am^ ^ 0. The numbering of massive neutrinos Vj is arbitrary. We 
will employ here the widely used convention of numbering of Vj which allows to associate #13 with 
the smallest mixing angle in the PMNS matrix U, and 612, Am^ > 0, and #23> Am^^j, with the 
parameters which drive, respectively, the solar [^] (z^ e ), and the dominant atmospheric (and 9^) 

2 Under the assumption of CPT invariance, which we will suppose to hold throughout this article, 812 and Am^i drive 
also the reactor v e oscillations at L ~ 180 km (see, e.g., [B]). 
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(and accelerator v^) oscillations. In this convention mi < 7712, < Am 2 i < |Am§i|, and, depending 
on sgn(Am3 1 ), we have either 7773 < m\ or 777,3 > 7772 (see further). In the case of m\ < 777,2 < 1^3 
(7773 < 777-1 < m 2)-, the neutrino mass squared difference Am^, as it follows from the data to be 
discussed below, is much smaller than lAm^^U, Am^ <C | Ant^/ggJ. This implies that in each of 
the two cases m x < m 2 < m 3 and m 3 < m x < m 2 we have |Am§! - Am 2 32 \ = Am 2 21 < |Am| 1)32 |. 
The angles #12 and 623 are sometimes called "solar" and "atmospheric" neutrino mixing angles, and 
are often denoted as #12 = 0© and 623 = a tm! while Am\ x and Am^jj are sometimes referred to 
as the "solar" and "atmospheric" neutrino mass squared differences and are sometimes denoted as |^] 
Am\ x = Ami, A77t| 1(32) = Am^ tm . 

The neutrino oscillation data, accumulated over many years, allowed to determine the parameters 
which drive the solar, reactor, atmospheric and accelerator neutrino oscillations, Am|i, ^12, lAmj^^J 
and 6*23, with a rather high precision (see, e.g., p2]). Furthermore, there were spectacular developments 
in the period since June 2011 in what concerns the angle #13 (see, e.g., [6]). They culminated in March 
of 2012 in a high precision determination of sin 2 20i3 in the Daya Bay experiment with reactor v e [18J: 



sin 2 26»i3 = 0.089 ± 0.010 ± 0.005 . (1.4) 

Subsequently]^ the RENO [H], Double Chooz, and T2K experiments [20] reported, respectively, 4.9<r, 
2.9o" and 3.2a evidences for a non-zero value of #13, compatible with the Day Bay result. 

A global analysis of the latest neutrino oscillation data presented at the Neutrino 2012 International 
Conference [T7], was performed in [21]. We give below the best fit values of Am|i, sin 2 612, lAm^^J, 
sin 2 623 and sin 2 #13, obtained in [2T] : 



Am\ x = 7.54 x 10" 5 eV 2 , |Am 2 1(32) | = 2.47 (2.46) x 10~ 3 eV 2 , (1.5) 
sin 2 0i2 = 0.307 , sin 2 23 = 0.39 , sin 2 6 13 = 0.0241 (0.0244) , (1.6) 

where the values (the values in brackets) correspond to mi < m2 < 7773 (777,3 < m i < mi)- The la 
uncertainties and the the 3a ranges of the neutrino oscillation parameters found in |21] are given in 
Table O 



A few comments are in order. We have Am| 1 /|Am 2 1 ^ 32 ^| = 0.031 <C 1, as was indicated earlier. 
The existing data do not allow to determine the sign of Am^^. As we will discuss further, the 
two possible signs correspond to two different basic types of neutrino mass spectrum. Maximal solar 
neutrino mixing, i.e. #12 = tt/4, is ruled out at more than 6a by the data. Correspondingly, one 



has cos 20i2 > 0.28 at 3a. The results quoted in eq. (1.6) imply that 023 — tt/4, 0i2 = 7r/5.4 and 
that 0i3 = 7r/20. Thus, the pattern of neutrino mixing is drastically different from the pattern of 
quark mixing. As we have noticed earlier, the neutrino oscillations experiments are sensitive only to 
neutrino mass squared differences Ara| = (m 2 — m 2 ), i 7^ j, and cannot give information on the 
absolute values of the neutrino masses, i.e., on the absolute neutrino mass scale. They are insensitive 
also to the nature - Dirac or Majorana, of massive neutrinos Vj and, correspondingly, to the Majorana 
CPV phases present in the PMNS matrix U [TTll2"2"]. 



3 The preceding part of the tex t of the present article follows closely parts of the text of the review article [B]. 

4 We have reported in eq. (1.4| the latest result of the Daya Bay experiment, published in the second article quoted 



in M, 

""' ' " ' ' ^"'31(32) 

is \Am\\ = |Arri§i - AmjL/21 (|Am 2 4 | = |AmI 2 + Am 2 21 /2\) 



5 Note that we have quoted the value of lAm^^)! in eq. (1.51, while the mass squared difference determined in [5T] 
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Table 1.1: The best-fit values and 3a allowed ranges of the 3- neutrino oscillation parameters, derived 
in [21] from a global fit of the current neutrino oscillation data. The values (values in brackets) 
correspond to mi < 777,2 < 7773 (7713 < mi < m 2 ). The definition of Am\ used is: Am\ = m\ — 
(m| + m 2 )/2. Thus, Am^ = Am^ — Am|i/2, if m\ < m 2 < 7713, and Am^ = Am^ + Am 2 i/2 for 

7773 < 771 1 < 777,2. 

Parameter best-fit (±ler) 3d 

Am| I 10 " 5 ey2 ] 
|Am|| [10" 3 eV 2 ] 
sin 2 6*12 
sin 2 6»23 
sin 2 6»i 3 



7.54^;^ 6.99 - 8.18 

2A3+m (2.42^;^) 2.19(2.17) - 2.62(2.61) 

0.307±g;gi| 0.259 - 0.359 

°- 386 iao2i (°- 392 io:oi) 0.331(0.335) - 0.637(0.663) 

0.0241±0.0025 (0.0244l°;°° 23 ) 0.0169(0.0171) - 0.0313(0.0315) 



After the successful measurement of #13, the determination of the absolute neutrino mass scale, of 
the type of the neutrino mass spectrum, of the nature - Dirac or Majorana, of massive neutrinos, as 
well as getting information about the status of CP violation in the lepton sector, are the most pressing 
and challenging problems and the highest priority goals of the research in the field of neutrino physics 
(see, e.g., [6j). 

As was already indicated above, the presently available data do not permit to determine the sign of 
Am^pj. In the case of 3-neutrino mixing, the two possible signs of Am 2 .-^^ correspond to two types 
of neutrino mass spectrum. In the widely used convention of numbering the neutrinos with definite 
mass employed by us, the two spectra read: 

i) spectrum with normal ordering (NO): mi < 7772 < 7773, Am 2 tm = Am^j > 0, Am?) = Am\i > 0, 

"72(3) = (™1 + A77l 2 1(31) )5 ; 

ii) spectrum with inverted ordering (IO): 7773 < 777-1 < m-2, Am 2 tm = A77732 < 0, At77q = AtttJi > 0, 

777 2 = (m 2 + A777 2 3 )l, 7771 = (777 2 + At77 2 3 - A777 2 i)i 

Depending on the value of the lightest neutrino mass, min(m./), the neutrino mass spectrum can be: 

a) Normal Hierarchical (NH): mi <^ m 2 < m 3 , m 2 = (Am|)^ = 8.68 x 10 -3 eV, 7773 = |Am 2 tm |5 ^ 
4.97 x 10~ 2 eV; or 

b) Inverted Hierarchical (IH): 7773 -C 7771 < 7772, with mi 2 — |Am 2 tm |a = 4.97 x 10 -2 eV; or 

c) Quasi-Degenerate (QD): mi = m 2 = 7773 = mo, m 2 3> |Am 2 tm |, mo *3 0.10 eV. 

The type of neutrino mass spectrum (hierarchy), i.e., the sign of AmL 32 j, can be determined i) 
using data from neutrino oscillation experiments at accelerators (NOz^A, T2K, etc.) (see, e.g., [23J), ii) 
in the experiments studying the oscillations of atmospheric neutrinos (see, e.g., (Sj)j as well as iii) in 
experiments with reactor antineutrinos |25j . The relatively large value of #13 is a favorable factor for 
the sgn(Am 2 1 ^ 32 - ) ) determination in these experiments. If neutrinos with definite mass are Majorana 
particles, information about the sgn(Am| 1 ^ 32 ^) can be obtained also by measuring the effective neutrino 
Majorana mass in neutrinoless double f3— decay experiments |26j . 

More specifically, in the cases i) and ii) the sgn(Am 31 ^ 32 -)) can be determined by studying oscilla- 
tions of neutrinos and antineutrinos, say, f„ v e and z/« -f> u e , in which matter effects are sufficiently 
large. This can be done in long base-line ^-oscillation experiments (see, e.g., [23]). For sin 2 2#i3 k; 0.05 
and sin 2 #23 ^ 0.50, information on sgn(Am 31 , 32 -,) might be obtained in atmospheric neutrino exper- 
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iments by investigating the effects of the sub-dominant transitions v^e) ~^ v e(n) an d ^W) — > ^e(n) °f 
atmospheric neutrinos which traverse the Earth (for a detailed discussion see, e.g., [24]). For v^ e ) 
(or ^V e )) crossing the Earth core, new type of resonance-like enhancement of the indicated transi- 
tions takes place due to the (Earth) mantle-core constructive interference effect (neutrino oscillation 
length resonance (NOLR)) [27j (see also [25]). As a consequence of this effect the corresponding v^ e ) 
(or ^( e )) transition probabilities can be maximal [29] (for the precise conditions of the mantle-core 
(NOLR) enhancement see [27J [29]) |^J For Am^^ > 0, the neutrino transitions v^e) ~^ v e(n) are 
enhanced, while for Am^pjj < the enhancement of antineutrino transitions 7^( e ) — > ^e(^) takes 
place [27] (see also [25J [291 [31] ) , which might allow to determine FJ sgn(Am| 1 ^ 32 ^). Determining the 
type of neutrino mass spectrum is crucial for understanding the origin of neutrino masses and mixing 
as well. 

All possible types of neutrino mass spectrum we have discussed above are compatible with the 
existing constraints on the absolute scale of neutrino masses mj. Information about the absolute 
neutrino mass scale can be obtained by measuring the spectrum of electrons near the end point in 
3 H /3-decay experiments [32] and from cosmological and astrophysical data (see, e.g., [33]). The most 
stringent upper bound on the v e mass was obtained in the Troitzk [35] experiment (see also |35j): 

m Pe < 2.05 eV at 95% C.L. (1.7) 

We have mp e = 7711,2,3 ^ 0.1 eV in the case of quasi-degenerate (QD) spectrum. The KATRIN 
experiment [35], which is under preparation, is planned to reach sensitivity of rrip e ~ 0.20 eV, i.e., 
it will probe the region of the QD spectrum [^] The Cosmic Microwave Background (CMB) data of 
the WMAP experiment, combined with supernovae data and data on galaxy clustering can be used 
to derive an upper limit on the sum of neutrinos masses (see, e.g., [33]). Depending on the model 
complexity and the input data used one obtains [37]: J2j m j & (0-3 — 1-3) eV, 95% C.L. Data on 
weak lensing of galaxies, combined with data from the WMAP and PLANCK experiments, may allow 
J2j mj to be determined with an uncertainty of cr(J2j m j) = (0.04 — 0.07) eV [38]. 

Thus, the data on the absolute scale of neutrino masses imply that neutrino masses are much 
smaller than the masses of the charged leptons and quarks. If we take as an indicative upper limit 
771 j ^ 0.5 eV, j = 1, 2, 3, we have 

— 3 — & 10~ 6 , I = e, fj,, r, q = d, s, b, u, c, t . (1-8) 
m i,q 

It is natural to suppose that the remarkable smallness of neutrino masses is related to the existence 
of a new fundamental mass scale in particle physics, and thus to new physics beyond that predicted 

6 We note that the Earth mantle-core (NOLR) enhancement of neutrino transitions differs [27] from the MSW one. It 
also differs |27II29] from the mechanisms of enhancement discussed, e.g., in the articles [3D]: the conditions of enhancement 
considered in [30] cannot be realised for the &Ve) ~~ * v e(n) or ^VO) — * ^e(m) transitions of the Earth core crossing neutrinos. 

7 The so-called "neutrino oscillograms of the Earth", showing in the case of Am^^) > 0, for instance, the dependence 
of, e.g., the VuM v e(p) transition probability on the nadir angle of the neutrino trajectory through the Earth for given 
values of the neutrino energy E (or Am^^-f/E) and of the relevant neutrino mixing angle, and discussed in detail in 
|31] and in the talks given by the authors of [31) . first appeared in the publications quoted in [28j . 

information on the type of neutrino mass spectrum can also be obtained in /3-decay experiments having a sensitivity 

to neutrino masses ~ Am^j^ | = 5 x f0 -2 eV [36] (i.e., by a factor of ~ 4 better sensitivity than that of the KATRIN 

experiment |35]V Reaching the indicated sensitivity in electromagnetic spectrometer /3-decay experiments of the type of 
KATRIN does not seem feasible at present. 
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by the Standard Model. A comprehensive theory of the neutrino masses and mixing should be able to 
explain the indicated enormous disparity between the neutrino masses and the masses of the charged 
leptons and quarks. 

At present no experimental information on the Dirac and Majorana CPV phases in the neutrino 
mixing matrix is available. Therefore the status of the CP symmetry in the lepton sector is unknown. 
The importance of getting information about the Dirac and Majorana CPV phases in the neutrino 
mixing matrix stems, in particular, from the possibility that these phases play a fundamental role in 
the generation of the observed baryon asymmetry of the Universe. More specifically, the CP violation 
necessary for the generation of the baryon asymmetry within the "flavoured" leptogenesis scenario 
|39[ 00] can be due exclusively to the Dirac and/or Majorana CPV phases in the PMNS matrix [41J, 
and thus can be directly related to the low energy CP-violation in the lepton sector. If the requisite CP 
violation is due to the Dirac phases 8, a necessary condition for a successful ("flavoured") leptogenesis 



is that sin #13 0.09 [41], which is comfortably compatible with the Daya Bay result, eq. (1.4). 

With ^ 0, the Dirac phase 8 can generate CP violating effects in neutrino oscillations [12] (see 
also [11, 43j), i.e., a difference between the probabilities of v\ — > vy and V[ — > vy oscillations in vacuum: 
P[yi — > vy) 7^ Pipi — > vy), I 7^ V = e,/x, r. The magnitude of the CP violating effects of interest is 
determined [3] by the rephasing invariant Jqp associated with the Dirac CPV phase 8 in U. It is 
analogous to the rephasing invariant associated with the Dirac CPV phase in the CKM quark mixing 



matrix |45j . In the "standard" parametrisation of the PMNS neutrino mixing matrix, eqs. (1.2) 



( 1.3 ), we have: 



J CP = Im (t^a U e2 U; 2 ) = i cos #13 sin 26 12 sin 2# 23 sin 26 13 sin 8. (1.9) 

Thus, given the fact that sin #12, sin #23 and sin #13 have been determined experimentally with a 
relatively good precision, the size of CP violation effects in neutrino oscillations depends essentially 
only on the magnitude of the currently unknown value of the Dirac phase 8. The current data imply 



\Jcp\ ^ 0.039, where we have used eq. (1.9) and the 3a ranges of sin 2 #12, sin 2 #23 and sin 2 #13 given 
in Table |1.1[ Data on the Dirac phase 8 will be obtained in the long baseline neutrino oscillation 
experiments T2K, NOz^A, etc. (see, e.g., refs. [6]). Testing the possibility of Dirac CP violation in the 
lepton sector is one of the major goals of the next generation of neutrino oscillation experiments (see, 
e.g., |23[ 146]). Measuring the magnitude of CP violation effects in neutrino oscillations is at present 
also the only known feasible method of determining the value of the phase 8 (see, e.g., |47|). 

If Vj are Majorana fermions, getting experimental information about the Majorana CPV phases 
in the neutrino mixing matrix U will be remarkably difficult |48|, |4"9"1 [50} [5T| [521 [53} [54] . As we will 
discuss further, the Majorana phases of the PMNS matrix play important role in the phenomenology 
of neutrinoless double beta ((/3/3)o u -) decay - the process whose existence is related to the Majorana 
nature of massive neutrinos [55]: (A, Z) — > (A, Z + 2) + e~ + e~. The phases 021,31 can affect 
significantly the predictions for the rates of the (LFV) decays \i — )• e + 7, r— >• /x + 7, etc. in a large 
class of supersymmetric theories incorporating the see-saw mechanism |56| . As was mentioned earlier, 
the Majorana phase(s) in the PMNS matrix can be the leptogenesis CP violating parameter(s) at the 
origin of the baryon asymmetry of the Universe. |4H 114] . 

Establishing whether the neutrinos with definite mass Vj are Dirac fermions possessing distinct 
antiparticles, or Majorana fermions, i.e., spin 1/2 particles that are identical with their antiparticles, 
is of fundamental importance for understanding the origin of i^-masses and mixing and the underlying 
symmetries of particle interactions. Let us recall that the neutrinos Vj with definite mass rrij will be 
Dirac fermions if particle interactions conserve some additive lepton number, e.g., the total lepton 
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charge L = L e + + L T . If no lepton charge is conserved, the neutrinos Vj will be Majorana fermions 
(see, e.g. [IS]). The massive neutrinos are predicted to be of Majorana nature by the see-saw mech- 
anism of neutrino mass generation [3] , which also provides an attractive explanation of the smallness 
of neutrino masses and, through the leptogenesis theory [2j, of the observed baryon asymmetry of the 
Universe. The observed patterns of neutrino mixing and of neutrino mass squared differences driving 
the solar and the dominant atmospheric neutrino oscillations, can be related to Majorana massive 
neutrinos and the existence of an approximate symmetry in the lepton sector corresponding to the 
conservation of the non-standard lepton charge L' = L e — — L T |57j . They can also be associated 
with the existence of approximate discrete symmetry (or symmetries) of the particle interactions (see, 
e.g., |58j). Determining the nature (Dirac or Majorana) of massive neutrinos is one of the fundamental 
and most challenging problems in the future studies of neutrino mixing |6j. 



1.2 The Nature of Massive Neutrinos 

1.2.1 Majorana versus Dirac Massive Neutrinos (Particles) 

The properties of Majorana particles (fields) are very different from those of Dirac particles (fields). 
A massive Majorana neutrino Xj ( or Majorana spin 1/2 particle) with mass rrij > can be described 
in local quantum field theory which is used to construct, e.g., the Standard Model, by 4-component 
complex spin 1/2 field Xj( x ) which satisfies the Dirac equation and the Majorana condition: 

C (x 3 ( x )) T = CkXj(x) , |&| 2 = 1, (1.10) 

where C is the charge conjugation matrix, C~ l j a C = — {^ a ) T (C T = — C, C _1 = and £j is, in 
general, an unphysical phase. The Majorana condition is invariant under proper Lorentz transforma- 
tions. It reduces by a factor of 2 the number of independent components in Xj( x )- 



The condition (1.10) is invariant with respect to U(l) global gauge transformations of the field 
Xj(x) carrying a U(l) charge Q, Xj{ x ) ~^ e%a Xji x ), om y if Q = 0. As a result, i) Xj cannot carry 
nonzero additive quantum numbers (lepton charge, etc.), and ii) the field Xj( x ) cannot "absorb" 
phases. Thus, Xj( x ) describes 2 spin states of a spin 1/2, absolutely neutral particle, which is identical 
with its antiparticle, Xj — Xj- As is well known, spin 1/2 Dirac particles can carry nonzero £7(1) 
charges: the charged leptons and quarks, for instance, carry nonzero electric charges. 

Owing to the fact that the Majorana (neutrino) fields cannot absorb phases, the neutrino mixing 
matrix U contains in the general case of n charged leptons and mixing of n massive Majorana neutrinos 
Uj = Xj, altogether 

(M) n ( n — 1) 

n CPV = 2 ' Majorana Uj , (1.1 1) 

CPV phases jTT]. In the case of mixing of n massive Dirac neutrinos, the number of CPV phases in 
U, as is well known, is 

"cpv = j . Dirac "i ■ ( 1 • 12 ) 

Thus, if Vj are Majorana particles, U contains the following number of additional Majorana CP 
violation phases: n^, p y = n^y — n|py = [n — 1). In the case of n charged leptons and n massive 
Majorana neutrinos, the PMNS matrix U can be cast in the form |11| 

U = VP, (1.13) 
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where the matrix V contains the (n — l)(n — 2)/2 Dirac CP violation phases, while P is a diagonal 
matrix with the additional (n — 1) Majorana CP violation phases ct2i, 0:31,..., ot n i, 



(■ a 21 ■ ( 
1, e l ^,e v 



■ 11 1 
2 



(1.14) 



As will discuss further, the Majorana phases will conserve CP if [15] ayi = irqj, qj = 0,1,2, j = 
2, 3, n. In this case exp(iaji) = ±1 and exp[z(aji — a^i)] = ±1 have a simple physical interpretation: 
these are the relative CP-parities of the Majorana neutrinos Uj and v\ and of Vj and u^, respectively. 

It follows from the preceding discussion that the mixing of massive Majorana neutrinos differs, in 
what concerns the number of CPV phases, from the mixing of massive Dirac neutrinos. For n = 3 
of interest, we have one Dirac and two Majorana CPV phases in U, which is consistent with the 
expression of U given in eq. (1.2). If n = 2, there is one Majorana CPV phase and no Dirac CPV 



phases in U. Correspondingly, in contrast to the Dirac case, there can exist CP violating effects even 
in the system of two mixed massive Majorana neutrinos (particles). 

The Majorana phases do not enter into the expressions of the probabilities of oscillations involving 
the flavour neutrinos and antineutrinos [1 1 1. 122j. v\ — > vy and V\ — > Vy. Indeed, the probability to find 
neutrino vy (antineutrino Vy) at time t if a neutrino vi (antineutrino V{) has been been produced at 
time to and it had traveled a distance L = t in vacuum, is given by (see, e.g., [TBI [6]): 



P(yi -)• vy 
P(pi ->vv) = 



j 



(1.15) 
(1.16) 



where Ej and pj are the energy and momentum of the neutrino Vj. It is easy to show, using the 
expression for U in eq. (1.13) that P[yi — > vy) and P(vi — > Vy) do not depend on the Majorana 
phases present in U since 



Vi e 



-i(Ej t-pj L) 



{vp)) l =Y,vv3^ i{Eit ' PiL) y 3 \ 



(1.17) 



The same result holds when the neutrino oscillations take place in matter |22j . 

If CP-invariance holds, Majorana neutrinos (particles) have definite CP-parity f]cp(xj) = i* 1 

Ucp Xj(x) Ucp = VCPiXj) 70 Xj(x P ), ricp(xj) = ipj = ±* , (1.18) 

where x = (xq,x), x p = (xq,— x) and Ucp is the unitary CP-transformation operator. In contrast, 
Dirac particles do not have a definite CP-parity — a Dirac field f(x) transforms as follows under the 
CP-symmetry operation: 



Ucpf(x)Ucp = r,f l0 C(f(x p )) T , 1 77/| 



1, 



(1.19) 



rjf being an unphysical phase factor. In the case of CP invariance, the CP-parities of massive Majorana 
fermions (neutrinos) can play important role in processes involving real of virtual Majorana particles 
(see, e.g., [El EDJ [BT]). 
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Using eqs. ( 1.19 ) and ( 1.18 ) and the transformation of the boson field under the CP-symmetry 
operation, 

U C p W a (x) U j CP = rj w K a (W a (xp))^ , \r] W \ 2 = 1 , « = -1 , «i,2,3 = +1 • (1-20) 

where r/w is an unphysical phase, one can derive the constraints on the neutrino mixing matrix U 
following from the requirement of CP-invariance of the leptonic CC weak interaction Lagrangian, eq. 
(1.1). In the case of massive Dirac neutrinos we obtain: rj^ViVwU^ = Uij, I = e,p,T, j = 1,2,3. 
Setting the product of unphysical phases rj*.r]irj]y = 1, one obtains the well known result: 



CP invariance 



Uij , I = e, p,,r , j = 1, 2, 3 , Dirac Vj 



(1.21) 



In the case of massive Majorana neutrinos we obtain using eqs. (1.10), (1.18) ( |1.19| ) and ( |1.20| ): 
S t j(ipj)r]ir]wUij = U^. It is convenient now to set £j = 1, rji = i and 7]w = 1. In this (commonly used 
by us) convention we get [To] : 



CP invariance : UZ = pj Uij , pj = +1 or (—1) , I = e, p,, r , j = 1, 2, 3 , Majorana Vj 



(1.22) 



Thus, in the convention used the elements of the PMNS matrix can be either real or purely imaginary 
if Vj are Majorana fermions. Applying the above conditions to, e.g., U e 2, U T 3 and U e ^ elements of the 



PMNS matrix (1.2) we obtain the CP conserving values of the phases 021, 031 and 5, respectively: 
«2i = kir, k = 0, 1, 2, 031 = k'ir, k! = 0, 1, 2, S = 0, ir, 2tt. 

One can obtain in a similar way the CP-invariance constraint on the matrix of neutrino Yukawa 
couplings, Xki, which plays a fundamental role in the leptogenesis scenario of baryon asymmetry 
generation, based on the (type I) see-saw mechanism of generation of neutrino masses [2"1 [3"1 ITT]: 



C Y (x) = - \ k i N kR {x) (x) ipi L {x) + h.c. , 
£m(z) = -\M k W k (x)N k {x). 



(1.23) 
(1.24) 



Here N k n{x) is the field of a heavy right-handed (RH) sterile Majorana neutrino with mass M k > 0, ipn 
denotes the Standard Model left-handed (LH) lepton doublet field of flavour / = e, p, r, rpf L = [yf L /J), 
and H is the Standard Model Higgs doublet field whose neutral component has a vacuum expectation 
value v = 174 GeV. The term Cy(x) + £^(x) includes all the necessary ingredients of the see-saw 
mechanism. Assuming the existence of two heavy Majorana neutrinos, i.e., taking k = 1,2 in eqs. 

to the Standard Model Lagrangian, we obtain 



(1.23) and (1.24), and adding the term C Y (x) + £ N 



M \ 



the minimal extension of the Standard Model in which the neutrinos have masses and mix and the 
leptogenesis can be realised. In the leptogenesis formalism it is often convenient to use the so called 
"orthogonal parametrisation " of the matrix of neutrino Yukawa couplings [59j: 



1 

v 



Am = - V M k R kj ^ffnTj {U ] )ji , 



(1.25) 



where R is, in general, a complex orthogonal matrix, R R T = R T R = 1. The CP violation necessary 
for the generation of the baryon asymmetry of the Universe is provided in the leptogenesis scenario 



of interest by the matrix of neutrino Yukawa couplings A (see, e.g., 0137]). It follows from eq. (1.25) 
that it can be provided either by the neutrino mixing matrix U, or by the matrix R, or else by both 
the matrices U and R. It is therefore important to derive the conditions under which A, R and U 



respect the CP symmetry. For the PMNS matrix U these conditions are given in eq. (1.22). For the 
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matrices A and R in the convention in which i) Nk(x) satisfy the Majorana condition with a phase 
equal to 1 (i.e., = 1), ii) rf = i and rj = 1, rf and r] H being the unphysical phase factors which 
appear in the CP-transformations of the LH lepton doublet and Higgs doublet fields [^J tPil(x) and 
H(x), respectively, they read [H] : 



A 



kl 



^kl Pk 



Pi" 



±1 , j = 1,2,3, I = e,/j,,r, 



R 



RkjPk Pj, 3, k = 1,2,3. 



(1.26) 
(1.27) 



where ip^ = rjcp(Nk) = ±i is the CP-parity of N^. Thus, in the case of CP invariance also the 



elements of A and R can be real or purely imaginary. Note that, as it follows from eqs. (1.22), (1.26) 



and (1.27), given which elements are real and which are purely imaginary of any two of the three 



matrices U, A and R, determines in the convention we are using and if CP invariance holds, which 
elements are real or purely imaginary in the third matrix. If, for instance, U e 2 is purely imaginary 
(p2 = —1) and Ai^ is real (p± = 1), then R\2 must be purely imaginary. Thus, in the example we are 
considering, a real R\2 would signal that the CP symmetry is broken |41j . 

The currents formed by Majorana fields have special properties, which make them also quite 
different from the currents formed by Dirac fields. In particular, it follows from the Majorana condition 
that the following currents of the Majorana field Xk( x ) are identically equal to zero (see, e.g., [To]): 

Xk{x)laXk{x) = 0, 

Xk{x)<y a i3 Xk(x) = 0, 
Xk{x)<y a pibXk{x) = 0. 



(1.28) 
(1.29) 
(1.30) 



Equations (1.28), (1.29) and (1.30) imply that Majorana fermions (neutrinos) cannot have nonzero 
U(l) charges and intrinsic magnetic and electric dipole moments, respectively. A Dirac spin 1/2 
particle can have non-trivial U(l) charges, as we have already discussed, and nonzero intrinsic magnetic 
moment (the electron and the muon, for instance, have it). If CP invariance does not hold, Dirac 



fermions can have also a nonzero electric dipole moments. Equations (1.29) and (1.30) imply also that 



the Majorana particles (neutrinos) cannot couple to a real photon. The axial current of a Majorana 
fermion, Xk{x) J a 75Xk(x) 7^ 0. Correspondingly, Xk(x) can have an effective coupling to a virtual 
photon via the anapole momentum term, which has the following form in momentum space: 



Qp) 7/3 75 



where q is the momentum of the virtual photon and Fa k \q' 2 ) is the anapole formfactor of Xk- The fact 
that the vector current of Xk is zero while the axial current is nonzero has important implications in 
the calculations of the relic density of the lightest and stable neutralino, which is a Majorana particle 
and the dark matter candidate in many supersymmetric (SUSY) extensions of the Standard Model 



(1.31) 



In certain cases (e.g., in theories with a keV mass Majorana neutrino (see, e.g., [9]), in the TeV 
scale type I see-saw model (see, e.g., [63]), in SUSY extensions of the Standard Model) one can have 
effective interactions involving two different massive Majorana fermions (neutrinos), say xi and X2- 



9 This convention is similar to, and consistent with, the convention about the unphysical phases we have used to derive 
the CP-invariance constraints on the elements of the PMNS matrix U. 



11 



Neutrino Physics 



The Nature of Massive Neutrinos 



We will consider two examples. The first is an effective interaction with the photon field, which can 
be written as: 

4$( x ) = Xi(*) °^(Pi2 - ^1275) X2(x) F a P(x) + h.c. , (1.32) 

where /ii2 and d\ 2 are, in general, complex constants, F a ^(x) = d a A^(x) — d^A a (x), A fl (x) being 
the 4- vector potential of the photon field. Using the Majorana conditions for Xi( x ) an d X2(x) in the 
convention in which the phases £1 = £2 = U it is not difficult to show that the constants p,\2 and 
d\2 enter into the expression for C^\x) in the form: (pi2 — /^* 2 ) = 2ilm(/xi 2 ) = /ii2, (^12 + d\ 2 ) = 
2Re(di 2 ) = d~i2, i-e., P12 is purely imaginary and J12 is real 10 In the case of CP invariance of 
the constants pu (A12) and d\2 (^12) should satisfy: 

CP invariance : P12 = - Pi P2 P12 , ^12 = + Pi P2 d\2 - (1.33) 

Thus, if pi = P2, i.e., if Xi{x) and X2(x) posses the same CP-parity, p\2 = and d\2 (and ^12) can be 
different from zero. If p\ = — P2, i-e., if Xi{ x ) and X2(x) have opposite CP-parities, d\2 = and p\2 
(and /i.12) can be different from zero. If CP invariance does not hold, we can have both p\2 7^ and 
di2 + (£12 + and d 12 + 0). 

As a second example we will consider effective interaction of xi an d X2 with a vector field (current), 
which for concreteness will be assumed to be the Z°-boson field of the Standard Model: 

d${x) = xi(x) la {v 12 - a 1275 ) X 2(x) Z a (x) + h.c. , (1.34) 

Here v\2 and a\2 are, in general, complex constants. Using the Majorana conditions for Xi( x ) and 
X2( x ) with £1 = £2 = U one can easily show that V12 has to be purely imaginary, while 012 has to be 
real 11 The requirement of CP invariance of £^(x), as can be shown, leads to (£1 = £2 = 1) : 

CP invariance : v± 2 = — pi P2 vu , ai 2 = + pi P2 ai2 ■ (1.35) 

Thus, we find, similarly to the case considered above, that if Xi(^) and X2(x) posses the same CP-parity 
(pi = P2), V12 = and a\2 can be different from zero; if Xi( x ) and X2(x) have opposite CP-parities 
(pi = — P2), a\2 = while v\2 can be different from zero. If CP invariance does not hold, we can have 
both v\2 7^ and a\2 / 0. 

These results have important implications, in particular, for the phenomenology of the heavy 
Majorana neutrinos in the TeV scale (type I) see-saw models, for the neutralino phenomenology in 
SUSY extensions of the Standard Model, in which the neutralinos are Majorana particles, and more 
specifically for the processes e~ + e + — > xi +X2-, X2 — > Xi + ^ + + ( m (X2) > rn{xi))-, I = P, T , where 
Xi and X2 are, for example, two neutralinos of, e.g., the minimal SUSY extension of the Standard 
Model (see, e.g., foOllol]). 

Finally, if ^/(x) is a Dirac field and we define the standard propagator of *&(x) as 

< 0\T(* a (x)*f,(y))\0 >= S^(x - y) , (1.36) 



10 In the case of Xi( x ) = Xz( x ) = x( x )> the current x{ x ) <T af){^i2 — di27s)X( :r ) nas to be hermitian, which implies 
that p,\2 should be real while di2 should be purely imaginary. Com bined with cons traints on /Z12 and di2 we have just 



obtained, this leads to /ii2 = di2 = 0, which is consistent with eqs. (1.291 and (1.301. 

n In the case of Xi( x ) = X^( x ) = x( x )i the hermiticity of the current x( x )la(vi2 — Qi2 / Yb)x( x ) implies that both V12 
and 012 have to be real. This, together with constraints on V12 and 012 just derived, leads to 1112 = 0, which is consistent 



with the result given in eq. ( 1.28 1 
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one has 

<0\T(* a (x)*f}(y))\0>=0 , < 0|T(# a (x)*f)(y))\0 >= . (1.37) 
In contrast, a Majorana neutrino field Xk{x) has, in addition to the standard propagator 

< 0\T( Xka (x)xkp(y))\Q >= Sj$(x - y) , (1.38) 
two non-trivial non-standard (Majorana) propagators 

< 0\T( X ka(x)xkp(y))\0 >= -&S%(x - y)C s ? , (1-39) 

< 0\T( Xka (x) X kp(y))\0 >= Cfe CjS^ix - y) . (1.40) 



This result implies that if Vj(x) in eq. ( 1.1 ) are massive Majorana neutrinos, (/3/3)o^-decay can proceed 
by exchange of virtual neutrinos Vj since < Q\T(yj a (x)i/jp(y))\Q >^ 0. The Majorana propagators play 
a crucial role in the calculation of the baryon asymmetry of the Universe in the leptogenesis scenario 
of the asymmetry generation (see, e.g., [3"1 14"?]). 

1.2.2 Generating Dirac and Majorana Massive Neutrinos 

The type of massive neutrinos in a given theory is determined by the type of the (effective) mass term 
C v m {x) neutrinos have, more precisely, by the symmetries of C v m {x) and of the total Lagrangian C{x) 
of the theory. A fermion mass term is any bi-linear in the fermion fields which is invariant under the 
proper Lorentz transformations. 

Massive Dirac neutrinos arise in theories in which the neutrino mass term conserves some additive 
quantum number that could be, e.g., the (total) lepton charge L = L e + L M + L T , which is conserved 
also by the total Lagrangian £(x) of the theory. A well known example is the Dirac mass term, which 
can arise in the minimally extended Standard Model to include three RH neutrino fields vm, I = e, fi, r, 
as SU{2)\ J singlets: 

L v D (x) = - Wr{x) M DVl u lL (x) + h.c. , (1.41) 

where Mq is a 3 x 3, in general complex, matrix. The term C v D {x) can be generated after the 
spontaneous breaking of the Standard Model gauge symmetry by an SU{2)l x U (1)y w invariant 
Yukawa coupling of the lepton doublet, Higgs doublet and the RH neutrino fields 



C Y (x) = -Y l Y l WR{x)Hl{x)tln L (x) + h.c, (1.42) 
M D = vY u . (1.43) 

If the nondiagonal elements of Mp are different from zero, Mdih ^ 0, I ^ I' = e, yu, r, the individual 
lepton charges Li, I = e, (jl, t, will not be conserved. Nevertheless, the total lepton charge L is conserved 
by C v D {x). As in the case of the charged lepton and quark mass matrices generated via the spontaneous 
electroweak symmetry breaking by Yukawa type terms in the SM Lagrangian, Mb is diagonalised by a 
bi-unitary transformation: Mq = U 1 ^ M D iag (U l £ p )\ where U l ^ p and U l ^ p are 3x3 unitary matrices. If 



the mass term in eq. (1.41 ) is written in the basis in which the charged lepton mass matrix is diagonal, 
U l ^ p coincides with the PMNS matrix, U^ p = t^PMNS- The neutrinos Uj with definite mass rrij > 
are Dirac particles: their fields Vjix) = (U^y^uiL^x) + (^jj P )]//^'_R(^) do not satisfy the Majorana 
condition, C {vj{x)) / S,jXji x )- Although the scheme we are considering is phenomenologically 
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viable 12 it does not provide an insight of why the neutrino masses are much smaller than the charged 
fermion masses. The only observable "new physics" is that related to the neutrino masses and mixing: 
apart from the neutrino masses and mixing themselves, this is the phenomenon of neutrino oscillations 



Indeed, given the fact that the lepton charges L\, I = e,\x,r, are not conserved, processes like 
fi + — > e + + 7 decay, pT — > e~ + e + + e~ decay, r~ — > e~ + 7 decay, etc. are allowed. However, the 
rates of these processes are suppressed by the factor EH \U V nU? A m 2 A I 'MlA 2 , V 4 I, Mw = 80 GeV 
being the W^-mass and I = [A, I' = e for the ^ — > & 
For instance, for the \i — > e + 7 decay branching ratio we have [61 



jVjufimj/M&l 2 , V + I, M w 
■ + 7 decay, etc., and are unobservably small. 



(see also |65j): 



BR(ji->e + -y) 



3a 
32n 



TJ ■ TJ* . I 

U V U M M 2 



W 



(2.5 - 3.9) x 10 



-55 



(1.44) 



where we have used the best fit values of the neutrino oscillation parameters given in eqs.(1.5) and 
( |1.6| ) and the two values correspond to 5 = tt and 0. The current experimental upper limit reads [66J: 
BR(fi + — > e + + 7) < 5.7 x 10 -13 . Thus, although the predicted branching ratio BR(fi + — > e + + 7) 7^ 0, 
its value is roughly by 43 orders of magnitude smaller than the sensitivity reached in the experiments 
searching for the \x — > e + 7 decay, which renders it unobservable in practice. 

As was emphasised already, massive Majorana neutrinos appear in theories with no conserved 
additive quantum number, and more specifically, in which the total lepton charge L is not conserved 
and changes by two units. In the absence of RH singlet neutrino fields in the theory, the flavour 
neutrinos and antineutrinos v\ and Pi, I = e,fj,,r, can have a mass term of the so-called Majorana 
type: 

£m0*0 = - I ^r(z) M in u lL (x) + h.c. , v c VR = C (PFl(x)) T , (1.45) 



2 " l ' Ry 



where M is a 3 x 3, in general complex, matrix. In the case when all elements of M are nonzero, 
Min 7^ 0, 1,1' = e,fj,,r, neither the individual lepton charges L; nor the total lepton charge L are 
conserved: L\ 7^ const., L ^ const. As it is possible to show, owing to the fact that un(x) are fermion 
(anti-commuting) fields, the matrix M has to be symmetric (see, e.g., |16|): M = M T . A complex 
symmetric matrix is diagonalised by the congruent transformation: 



M diag _ tj^mU , U - unitary . 



(1.46) 



where U is a 3 x 3 unitary matrix. If C v M {x) is written in the basis in which the charged lepton mass 
matrix is diagonal, U coincides with the PMNS matrix: U = C^pmns- The fields of neutrinos Uj with 
definite mass nij are expressed in terms of and uf R : 



Vj(x) 



-u j (x)m j u j {x) , 



u} lVlL {x) + Ulv c lR = C (UJ(x)) t , j = 1,2,3. 



(1.47) 
(1.48) 



They satisfy the Majorana condition with £,■ = !, as eq. (1.48) shows 



The Majorana mass term (1.45) for the LH flavour neutrino fields vil can be generated 

i) effectively after the electroweak symmetry (EWS) breaking in the type I see-saw models [I], 

ii) effectively after the EWS breaking in the type III see-saw models [6 

12 It does not contain a candidate for a dark matter particle. 
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hi) directly as a result of the EWS breaking by an SU(2) l triplet Higgs held which carries two units of 
the weak hypercharge Y\y and couples in an SU(2)l x Uy w invariant manner to two lepton doublets 

(the Higgs Triplet Model (HTM) sometimes called also "type II see-saw model"), 
iv) as a one loop correction to a Lagrangian which does not contain a neutrino mass term [69J (see 
also [7U]), 

iv) as a two loop correction in a theory where the neutrino masses are zero at tree and one loop levels 
[7T] (see also [70]), 

v) as a three loop correction in a theory in which the neutrino masses are zero at tree, one loop and 
two loop levels [70] . 

In all three types of see-saw models, for instance, the neutrino masses can be generated at the EWS 
breaking scale and in this case the models predict rich beyond the Standard Model physics at the TeV 
scale, some of which can be probed at the LHC (see, e.g., [72] and further). We will consider briefly 
below the neutrino mass generation in the type I see-saw and the Higgs triplet models. 



In a theory in which the SU{2)i J singlet RH neutrino fields vm, I = e, /i, r, are present , the most 



general neutrino mass Lagrangian contains the Dirac mass term (1.41), the Majorana mass term for 



the LH flavour neutrino fields (1.45) and a Majorana mass term for the RH neutrino fields vir(x) |74j : 



£-d+m( x ) = ~ vvr(x) Mni'i vil(x)- - vf, R (x) M LV \ vil(x)- - vf, L (x) M RVi v m (x) + h.c. , (1.49) 

where uf, L = C (ujT R (x)) T and Mp, Ml and M R are 3 x 3, in general complex, matrices. By a simple 
rearrangement of the neutrino fields this mass term can be cast in the form of a Majorana mass term 
which is then diagonalised with the help of the congruent transformation |16j . In this case there are 
six Majorana mass eigenstate neutrinos, i.e., the flavour neutrino fields vil(x) are linear combinations 



of the fields of six Majorana neutrinos with definite mass. The neutrino mixing matrix in eq. (1.1) is 
a 3 x 6 block of a 6 x 6 unitary matrix. 

The Dirac-Majorana mass term is at the basis of the type I see-saw mechanism of generation of 
the neutrino masses and appears in many Grand Unified Theories (GUTs) (see, e.g., |16| for further 
details). In the see-saw models, some of the six massive Majorana neutrinos typically are too heavy to 
be produced in the weak processes in which the initial states of the flavour neutrinos and antineutrinos 
vi and used in the neutrino oscillation experiments, are being formed. As a consequence, the states 
of vi and vg will be coherent superpositions only of the states of the light massive neutrinos fj, and 
the elements of the neutrino mixing matrix C/pmnSj which are determined in experiments studying 
the oscillations of vg and will exhibit deviations from unitarity. These deviations can be relatively 
large and can have observable effects in the TeV scale see-saw models, in which the heavy Majorana 
neutrinos have masses in the ~ (100 — 1000) GeV range (see, e.g., [75]). 

If after the diagonalisation of C U D+M {x) more than three neutrinos will turn out to be light, i.e., to 
have masses ~ 1 eV or smaller, active-sterile neutrino oscillations can take place (see, e.g., (9[[T6j]): a 
LH (RH) flavour neutrino vil (antineutrino vi R ) can undergo transitions into LH sterile antineutrino(s) 
vi'L = v s vl (RH sterile neutrino(s) vy R = uf, R ). As a consequence of this type of oscillations, one 
would observe a "disappearance" of, e.g., v e and/or z/« (p e and/or z/u) on the way from the source to 
the detector. 

We would like to discuss next the implications of CP invariance for the neutrino Majorana mass 



matrix, eq. (1.45). In the convention we have used to derive eqs. (1.26) and (1.27), in which the 



13 We consider in the present article the case of three RH sterile neutrinos, but schemes with less than 3 and more than 
3 sterile neutrinos are also discussed in the literature, see, e.g., [9]. 
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unphysical phase factor in the CP transformation of the lepton doublet field ^il(x), and thus of 
vil{x), rf = i, the requirement of CP invariance leads to the reality condition for M: 



CP 



invariance 



M* = M. 



(1.50) 



Thus, M is real and symmetric and therefore is diagonalised by an orthogonal transformation, i.e., if 



CP invariance holds, the matrix U in eq. (1.46) is an orthogonal matrix. The nonzero eigenvalues of 
a real symmetric matrix can be positive or negative^] Consequently, M d%ag in eq. (1.46) in general 
has the form: 



M dia9 = (mi,m' 2 ,m£), m' 



m 



> , pj = ±l. 



ri.5i 



Let us denote the neutrino field which has a mass m'- ^ by v'Ax). According to eq. (1.48), the 



field v'jix) satisfies the Majorana condition: v',-{x) = C (z/(x)) T . One can work with the fields v'j(x) 
remembering that some of them have a negative mass. It is not difficult to show that the CP-parity 



of the fields v'Ax) is r)cp{v'j] 



h J 



1,2,3. The physical meaning of the signs of the masses m'- ^ 



of the Majorana neutrinos becomes clear if we change to a "basis" of neutrino fields Vj(x) which have 
positive masses rrij > 0. This can be done, e.g., by introducing the fields [TE] : 



v 'ji. x ) = (~75)2 (1 pj) Vj(x) : Vj(x) = Uj{x) if pj = 1 ; Vj(x) = -75 Vj(x) if 



Pi 



1 . 



(1.52) 



As it is not difficult to show, if v'Ax) has a mass m'j < 0, CP-parity r)cp(Vj) = i and satisfies 
the Majorana condition C (fj(x)) T = pjis'j(x), the field Vj(x) posses a mass rrij > 0, CP-parity 
Vcpi^j) = ipj, and satisfy the Majorana condition C (l7j(x)) T = pji/j(x): 

Vj : rrij > 0, r]cp(vj) = ipj , C (^"(x)) T = pjVj(x) . (1.53) 

Thus, in the case of CP invariance, the signs of the nonzero eigenvalues of the neutrino Majorana 
mass matrix determine the CP-parities of the corresponding positive mass Majorana (mass eigenstate) 
neutrinos 15 



1.2.3 A Brief Historical Detour 

It is interesting to note that B. Pontecorvo in his second seminal article on neutrino oscillations [7], 
which was published in 1958 when only one type of neutrino and antineutrino was known, assumed 
that the state of the neutrino v, emitted in weak interaction processes is a linear superposition of 
the states of two Majorana neutrinos v^ 1 and vtf which have different masses, m\ ^ 1712, opposite 
CP-parities, ??cp( i/ i V/ ) = ~ t ]cp{ ii 2 1 ) an d are maximally mixed, while the state of the corresponding 
antineutrino V is just the orthogonal superposition of the states of and v^: 

\u> = J , (1.54) 

\u M > - \ u M > 

\v> = J ■ (1.55) 

14 The absolute value of the difference between the number of positive and number of negative eigenvalues of a real 
symmetric matrix A is an invariant of the matrix with respect to transformations A' — PAP T , where P is a real matrix 
which has an inverse. 

15 For further discussion of the properties of massive Majorana neutrinos (fermions) and their couplings see, e.g., [16] . 
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Thus, the oscillations are between the neutrino v and and the antineutrino V, in full analogy with 
the K° — K° oscillations. From contemporary point of view, B. Pontecorvo proposed active - sterile 
neutrino oscillations with maximal mixing and massive Majorana neutrinos. To our knowledge, the 
1958 article quoted in [7] was also the first in which fermion mixing in the weak interaction Lagrangian 
was introduced. 

The article of Z. Maki, M. Nakagawa and S. Sakata [8] was inspired, in part, by the discovery of the 
second type of neutrino - the muon neutrino, in 1962 at Brookhaven. These authors considered a com- 
posite model of elementary particles in which the electron and muon neutrino states are superpositions 
of the states of composite Dirac neutrinos and which have different masses, mf 7^ m®: 

\v e > = \ui > cos 9 C + \v§ > sin 9 C , (1.56) 
> = - |z/f > sin 9 C + \i>2 > cos 9 C , (1-57) 

where 9 C is the neutrino mixing angle. The model proposed in [8] has lepton-hadron symmetry built 
in and as consequence of this symmetry the neutrino mixing angle coincides with what we call today 



the Cabibbo angle 9 C = 0.22. The authors of [8] discuss the possibility of — v e oscillations, which 



they called "virtual transmutations" . 

In an article [75] by Y. Katayama, K. Matsumoto, S. Tanaka and E, Yamada, published in 1962 
somewhat earlier than [8], the authors also introduce two-neutrino mixing. However, this is done purely 
for model construction purposes and does not have any physical consequences since the neutrinos in 
the model constructed in |78j are massless particles. 

In 1967 B. Pontecorvo independently considered the possibility of v e z/„ oscillations in the article 
PP, in which the notion of a "sterile" or "inert" neutrino was introduced. Later in 1969, V. Gribov and 
B. Pontecorvo [79] introduced for the first time a Majorana mass term for the LH flavour neutrinos 
v e and i/„, the diagonalisation of which lead to two Majorana neutrinos v^2 with definite but different 
masses mi^, m\ 7^ m,2, and two- neutrino mixing with an arbitrary mixing angle 9: 

\u e > = \ v f > cos0 + \vf > sin9, (1.58) 
1 1/,, > = - \vi > sine + \vf > cos 9. (1.59) 

This was the first modern treatment of the problem of neutrino mixing which anticipated the way 
this problem is addressed in gauge theories of electroweak interactions and in Grand Unified Theories 
(GUT's). In the same article for the first time the analytic expression for the probability of u e o 
oscillations was also derived. 

1.2.4 Models of Neutrino Mass Generation: Two Examples 

Type I See-Saw Model. A natural explanation of the smallness of neutrino masses is provided by 
the type I see-saw mechanism of neutrino mass generation [3]. An integral part of this rather simple 
mechanism are the RH neutrinos vm (RH neutrino fields vm(x)). The latter are assumed to possess a 
Majorana mass term £^j(x) as well as Yukawa type coupling £y(x) with the Standard Model lepton 



and Higgs doublets, i/jil(x) and H(x), given in eq. (1.42). In the basis in which the Majorana mass 
matrix of RH neutrinos is diagonal, we have: 

C YM (x) = C Y (x) + C^(x) = - (a h AWx) H\x) ih L (x) + h.c.) - X - M k W k (x) N k {x) , (1.60) 



16 The article by Z. Maki, M. Nakagawa and S. Sakata [8] appeared before the article by N. Cabibbo [76] in which the 
"Cabibbo angle" 9 C was introduced and the hadron phenomenology related to this angle was discussed, but after the 
article by M. Gell-Mann and M. Levy [77] in which 9 C was also introduced (by the way, in a footnote). 
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where we have combined the expressions given in eqs. (1.23) and (1.24). When the electroweak 



symmetry is broken spontaneously, the neutrino Yukawa coupling generates a Dirac mass term: 
m M Nkn(x) vil(x) +h.c, with mP = vX, v = 174 GeV being the Higgs doublet v.e.v. In the case when 
the elements of mP are much smaller than M k , \mP, \ <C i, k = 1, 2, 3, I = e,^,r, the interplay 
between the Dirac mass term and the mass term of the heavy (RH) Majorana neutrinos Nj- generates 
an effective Majorana mass (term) for the LH flavour neutrinos [3]: 

M n - - {m D )J, k Mfr 1 m° = -v 2 (X)J! k M k l X kl . (1.61) 

In grand unified theories, m D is typically of the order of the charged fermion masses. In SO(W) 
theories, for instance, mP coincides with the up-quark mass matrix. Taking indicatively M ~ 0.1 
eV, m D ~ 100 GeV, one obtains M k ~ Mjq ~ 10 14 GeV, which is close to the scale of unification of 
the electroweak and strong interactions, Mqut — 2 x 10 16 GeV. In GUT theories with RH neutrinos 
one finds that indeed the heavy Majorana neutrinos N k naturally obtain masses which are by few to 
several orders of magnitude smaller than Mqut (see, e.g., the second and third articles quoted in ref. 
[1]). Thus, the enormous disparity between the neutrino and charged fermion masses is explained in 
this approach by the huge difference between effectively the electroweak symmetry breaking scale and 
Mqut- 

An additional attractive feature of the see-saw scenario under discussion is that the generation and 
smallness of neutrino masses is related via the leptogenesis mechanism [2] (see also, e.g., [3| [39 1 l4"0l I4"T] ) 



to the generation of the baryon asymmetry of the Universe. Indeed, the Yukawa coupling in ( 1.60), in 
general, is not CP conserving. Due to this CP-nonconserving coupling the heavy Majorana neutrinos 
undergo, e.g., the decays Nj — > l + + H^~\ Nj — > l~ + H^ + \ which have different rates: T(Nj — > 
l + + H^P) / T{Nj -»• r + # (+) )- When these decays occur in the Early Universe at temperatures 
somewhat below the mass of, say, Aq, so that the latter are out of equilibrium with the rest of the 
particles present at that epoch, CP violating asymmetries in the individual lepton charges L\ and 
in the total lepton charge L of the Universe are generated. These lepton asymmetries are converted 
into a baryon asymmetry by (B — L) conserving, but (B + L) violating, sphaleron processes, which 
exist in the Standard Model and are effective at temperatures T ~ (100 - 10 12 ) GeV [80]. If the 
heavy neutrinos Nj have hierarchical spectrum, M\ <C M2 <C M3, the observed baryon asymmetry 
can be reproduced provided the mass of the lightest one satisfies Mi <; 10 9 GeV [HI] ^\ Thus, in 
this scenario, the neutrino masses and mixing and the baryon asymmetry have the same origin - the 
neutrino Yukawa couplings and the existence of (at least two) heavy Majorana neutrinos. Moreover, 
quantitative studies based on advances in leptogenesis theory |39| I40j. in which the importance of the 
flavour effects in the generation of the baryon asymmetry was understood, have shown that the Dirac 
and/or Majorana phases in the neutrino mixing matrix U can provide the CP violation, necessary in 
leptogenesis for the generation of the observed baryon asymmetry of the Universe |41j . This implies, 
in particular, that if the CP symmetry is established not to hold in the lepton sector due to the PMNS 
matrix U, at least some fraction (if not all) of the observed baryon asymmetry might be due to the 
Dirac and/or Majorana CP violation present in the neutrino mixing. 

In the see-saw scenario considered, the scale at which the new physics manifests itself, which is 
set by the scale of masses of the RH neutrinos, can, in principle, have an arbitrary large value, up to 
the GUT scale of 2 x 10 16 GeV and even beyond, up to the Planck mass. An interesting possibility, 
which can also be theoretically well motivated (see, e.g., |83j), is to have the new physics at the TeV 
scale, i.e., M k ~ (100 — 1000) GeV. Low scale see-saw scenarios usually predict a rich phenomenology 

17 In specific type I see-saw models this bound can be lower by a few orders of magnitude, see, e.g., |82| . 
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at the TeV scale and are constrained by different sets of data, such as, e.g., the data on neutrino 
oscillations, from EW precision tests and on the lepton flavour violating (LFV) processes [i —> cy, 
/j, — > 3e, pT — e~ conversion in nuclei. In the case of the TeV scale type I see-saw scenario of interest, 
the flavour structure of the couplings of the heavy Majorana neutrinos Nk to the charged leptons and 
the bosons, and to the LH flavour neutrinos un and the Z° boson, are essentially determined 
by the requirement of reproducing the data on the neutrino oscillation parameters [63J. All present 
experimental constraints on this scenario still allow i) for the predicted rates of the fj, — )■ e + 7 decay, 
fi — > 3e decay and fj, — e conversion in the nuclei to be [S3] within the sensitivity range of the currently 
running MEG experiment on /x — > e+7 decay [?] planned to probe values of BR(n + — > e + +7) ^ 10~ 13 , 
and of the future planned experiments on \i — > 3e decay and \x — e conversion [85], [86, 87, 88, 89], h) for 
an enhancement of the rate of neutrinoless double beta ((/3/3)o v -) decay [63], which thus can be in the 
range of sensitivity of the (/3/3)oi/-decay experiments which are taking data or are under preparation 
(see, e.g., |90j ) even when the light Majorana neutrinos possess a normal hierarchical mass spectrum 
(see further), and iii) for the possibility of an exotic Higgs decay channel into a light neutrino and a 
heavy Majorana neutrino with a sizable branching ratio, which can lead to observable effects at the 
LHC |72j (for further details concerning the low energy phenomenology of the TeV scale type I see-saw 
model see, e.g., [631 183] fM]V 

Let us add that the role of the experiments searching for lepton flavour violation to test and possibly 
constrain low scale see-saw models, and more generally, extensions of the Standard Model predicting 
"new" (lepton flavour violating) physics at the TeV scale, will be significantly strengthened in the 
next years. Searches for fj, — e conversion at the planned COMET experiment at KEK [86] and Mu2e 
experiment at Fermilab [57] aim to reach sensitivity to conversion rates CR(/xAl — > eAl) ~ 10 -16 , 
while, in the longer run, the PRISM/PRIME experiment in KEK [88] and the project-X experiment 
in Fermilab [89J are being designed to probe values of the \i — e conversion rate on Ti, which are by 
2 orders of magnitude smaller, CR(/iTi — > eTi) ~ 10 -18 [88| . The current upper limit on the [i — e 
conversion rate is CR(//A1 — > eAl) < 4.3 x 1CP 12 [91]. There are also plans to perform a new search 
for the /i + — > e + e~e + decay |85J, which will probe values of the corresponding branching ratio down 
to BR(^ + — > e + e~e + ) « 10~ 15 , i.e., by 3 orders of magnitude smaller than the best current upper 
limit [92] . Furthermore, searches for tau lepton flavour violation at superB factories aim to reach a 
sensitivity to BR(r — > (/j,,e)"f) « 10~ 9 (see, e.g., |93j). 

The Higgs Triplet Model (HTM). In its minimal formulation this model includes one addi- 
tional SU{2)l triplet Higgs field A, which has weak hypercharge Yiy = 2 [68J: 



The Lagrangian of the Higgs triplet model which is sometimes called also the "type II see-saw model" , 



where ip C £L = (— f^C -1 — £^C _1 ), C being the charge conjugation matrix, H is the SM Higgs 
doublet and /Ua is a real parameter characterising the soft explicit breaking of the total lepton charge 
conservation. We will discuss briefly the low energy version of HTM, where the new physics scale 

18 We do not give here, for simplicity, all the quadratic and quartic terms present in the scalar potential (see, e.g., |95]~). 




(1.62) 



reads 



£htm 



Ml Tr (At A) - (hu, i> c tL ir 2 Ai/je L + fi A H ] A f ir 2 H* + h.c.) , (1.63) 
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Ma associated with the mass of A takes values 100 GeV & Ma ^ 1 TeV, which, in principle, can be 
probed by LHC (see [HUES] and references quoted therein). 

The flavour structure of the Yukawa coupling matrix h and the size of the lepton charge soft 
breaking parameter //a are related to the light neutrino Majorana mass matrix M v , which is generated 
when the neutral component of A develops a "small" vev va cx pa • Indeed, setting A = va and 



H T = (v 0) T with v ~ 174 GeV, from Lagrangian (1.63) one obtains: 



M v u , ~ 2h U 'V A . (1.64) 

The matrix of Yukawa couplings hw is directly related to the PMNS neutrino mixing matrix C/pmns = 
U, which is unitary in this case: 

h u > = ^ ( U * diag(mi,m 2 ,m 3 ) ^ f ) a , , mj > . (1.65) 

An upper limit on va can be obtained from considering its effect on the parameter p = M^ /M| cos 2 9yy. 
In the SM, p = 1 at tree-level, while in the HTM one has 

1 + 2x 2 

p=l + 5p= - - ^ 2 , x = v A /v. (1.66) 

The measurement p « 1 leads to the bound va/v ^ 0.03, or v a < 5 GeV (see, e.g., [95]). 

For Ma ~ (100 — 1000) GeV, the model predicts a plethora of beyond the SM physics phenomena 
(see, e.g., (95j EH EZ]), most of which can be probed at the LHC and in the experiments on charged 
lepton flavour violation, if the Higgs triplet vacuum expectation value va is relatively small, roughly 
va ~ (1 — 100) eV. As can be shown (see, e.g., [95]), the parameters va and pa are related: for 
Ma ~ v = 174 GeV we have va = PA, while if M A >> v 2 , then v A = pav 2 /(2M a ). Thus, a 
relatively small value of va in the TeV scale HTM implies that pa has also to be small, and vice verse. 
A nonzero but relatively small value of pa can be generated, e.g., at higher orders in perturbation 
theory [98] . The smallness of the neutrino masses is therefore related to the smallness of the vacuum 
expectation value va, which in turn is related to the smallness of the parameter pa- 

Under the conditions specified above one can have testable predictions of the model in low energy 
experiments, and in particular, in the ongoing MEG and the planned future experiments on the lepton 
flavour violating processes p — > ej, p — > 3e and p + M — > e + M (see, e.g., [M]). The HTM has also 
an extended Higgs sector including neutral, singly charged and doubly charged Higgs particles. The 
physical singly-charged Higgs scalar field (particle) practically coincides with the triplet scalar field 
A + , the admixture of the doublet charged scalar field being suppressed by the factor va/v. The singly- 
and doubly- charged Higgs scalars A + and A ++ have, in general, different masses [98 1: rn-A+ 7^ m,A++- 
Both cases toa+ > wa++ an d niA+ < ?71a++ are possible. The TeV scale HTM predicts the existence 
of rich new physics at LHC as well, associated with the presence of the singly and doubly charged 
Higgs particles A + and A ++ in the theory (see, e.g., [95, 97]). 



1.3 Determining the Nature of Massive Neutrinos 

The Majorana nature of massive neutrinos typically manifests itself in the existence of processes 
in which the total lepton charge L changes by two units: K + — > ir~ + p + + p + , pT + (A, Z) — > 
p + + (A, Z — 2), etc. Extensive studies have shown that the only feasible experiments having the 
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potential of establishing the Majorana nature of massive neutrinos at present are the (/3/3)oi/-decay 
experiments searching for the process (A, Z) — > (A, Z + 2) + e~ + e~ (for reviews see, e.g., [To] [90 ] [99] ) . 
The observation of (/3/3)ojy-decay and the measurement of the corresponding half-life with sufficient 
accuracy, would not only be a proof that the total lepton charge is not conserved, but might provide 
also information i) on the type of neutrino mass spectrum J2H] , and ii) on the absolute scale of neutrino 
masses (see, e.g. [4"9]). 

The observation of (/3/3)oi/-decay and the measurement of the corresponding half-life with sufficient 
accuracy, combined with data on the absolute neutrino mass scale might provide also information on 
the Majorana phases in U |100] B£>, 49, 50, 5JJ- If the neutrino mass spectrum is inverted hierarchical or 
quasi-degenerate, for instance, it would be possible to get information about the phase 021- However, 
establishing even in this case that 021 has a CP violating value would be a remarkably challenging 
problem [51] (see also [52]). Determining experimentally the values of both the Majorana phases 
021 an d 031 is an exceptionally difficult problem. It requires the knowledge of the type of neutrino 
mass spectrum and high precision determination of both the absolute neutrino mass scale and of the 
(/3/3)oi/-decay effective Majorana mass, |( m )| (see, e.g., [4"8 l [50 ] [5T]). 



1.3.1 Majorana Neutrinos and (/3(3) 0v — Decay- 
Under the assumptions of 2>-v mixing, for which we have compelling evidence, of massive neutrinos 
Uj being Majorana particles, and of (/3/3)oi/-decay generated only by the (V-A) charged current weak 
interaction via the exchange of the three Majorana neutrinos Vj having masses rrij ^ few MeV, the 
(/3/3)oi/-decay amplitude of interest has the form (see, e.g. [TBI SSI HS3 ) : A(f3f3)o u = < m > M, where 
M is the corresponding nuclear matrix element (NME) which does not depend on the neutrino mixing 
parameters, and 

|(m)| = \nn\Ua\ 2 + m 2 \U e2 \ 2 e ia21 +m 3 \U e3 \ 2 e i ^ 1 - 2S ^\ , (1.67) 

is the effective Majorana mass in (/?/3)ot/-decay, \U e i\ = C12C13, \U e 2\ = S12C13, \U e3 \ = S13. In the case 
of CP-invariance one has 25 = or 2ir and [15} [13] , 

r?2i = e* Q21 = ±1 , r/31 = e ia ^ = ±1 , (1.68) 

r?2i(3i) being the relative CP-parity of the Majorana neutrinos ^2(3) and v\. 

It proves convenient to express [101] the three neutrino masses in terms of Am^ and Am 2 -^^^, 
measured in neutrino oscillation experiments, and the absolute neutrino mass scale determined by 
min(mj) [^J In both cases of neutrino mass spectrum with normal and inverted ordering one has 
(in the convention we use): Am^ > 0, rri2 = (m 2 + Am^)^ For normal ordering, Ami > and 
TO3 = (m 2 + Am 2 ^! , while if the spectrum is with inverted ordering, min(mj) = m 3 , Am 2 2 < and 
m\ = (m 2 + A77123 — A m 2i) 2 • Thus, given Am^, Am^^, 6*12 and #13, |( m )\ depends on min(rrij), 
Majorana phases 021, 031 and the type of v-mass spectrum. 

The problem of obtaining the allowed values of |( m }| given the constraints on the parameters 
following from z/-oscillation data, and more generally of the physics potential of (/3/3)o^- decay ex- 



periments, was first studied in [101J and subsequently in a large number of papers (see, e.g. 



[4"5T EJJ 1102} 11031 1104] ) . The results of this analysis are illustrated in Fig. 1 . The main features of the 



19 For a detailed discussion of the relevant formalism see, e.g. (161 1481 [99] . 
20 Extensive list of references on the subject is given in [99] . 
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predictions for | ( m ) | in the cases of the NH, IH and QD spectra are summarised below, 
i) NH spectrum: 

|< m )| = |(Am2i)2 s\ 2 + (AmlOi s 2 3 e^-^^ . (1.69) 
Using the 3<r allowed ranges of the relevant neutrino oscillation parameters we get: 

4.7 x 1(T 4 eV £ |( m )\ & 4.8 x l(T 3 eV , NH . (1.70) 



ii) IH spectrum: 



Numerically one finds: 



K ™ }| = (|AmI a |)3 (1 - sin 2 2# 21 sin 2 ^ , (1.71) 
(|Am| 2 |)3 cos2(9i2 £ |( m )| £ (|Am| 2 |)5 . (1.72) 



0.014 eV £ |( m )| £ 0.050 eV , IH , (1.73) 



the upper and the lower bounds corresponding to the CP-conserving values of a 2 \ = 0; n. 
hi) QD spectrum: 

|( m >| ^m (l-sin 2 26> 12 sin 2 ^)5 , (1.74) 
m |( m )| £ m cos20i 2 £ 0.028 eV , (1.75) 



with m £ 0.1 eV, m < 2.05 eV [M] (see also [35]), or m & (0.3 - 1.3) eV [37] (see eq. (JL7J) and 
the discussion following after it). 



For the IH (QD) spectrum we have also [100, 

I 

2 ' V rh' 2 I sin^20i 2 



sin 2 (°^) - 1 - , m 2 s | Am 2 2 | (m 2 } f m (QD) _ (L?6) 



Thus, a measurement of |( m )| and mo (lAm^l) for QD (IH) spectrum can allow to determine a 2 \. 

The experimental searches for (/3/3)o;,-decay have a long history (see, e.g., [105j). The most stringent 
upper limits on |( m )| were set by the IGEX [TO6], CUORICINO [TOT] . NEM03 [IDS] and EXO-200 
[T09] experiments with 76 Ge, 130 Te, 100 Mo and 136 Xe, respectively The IGEX collaboration has 
obtained for the half-life of 76 Ge, 1% > 1.57 x 10 25 yr (90% C.L.), from which the limit |( m )\ < 
(0.33 — 1.35) eV was derived [106] . Using the recent more advanced calculations of the corresponding 
nuclear matrix elements (including the relevant uncertainties) jllOj one finds: |( m )| < (0.22 — 
0.35) eV. The NEM03 and CUORICINO experiments, designed to reach a sensitivity to |( m )| ~ 
(0.2 - 0.3) eV, set the limits: |( m }| < (0.61 - 1.26) eV [JUS] and |( m )| < (0.19 - 0.68) eV [107] 
(90% C.L.), where estimated uncertainties in the NME are accounted for. The two upper limits were 
derived from the experimental lower limits on the half-lives of 100 Mo and 130 Te, Ty 2 > 5.8 x 10 23 yr 

^ q n \, i n24 Tr „ (c\c\o/.n t \ Ii n>7i w;+v. +r,^ wmt?„ 

1/2 



C.L.) [TUB] and Ty 2 > 3.0 x 10 24 yr (90%C.L.) [107J. With the NMEs and their uncertainties 



lr The NEM03 collaboration has searched for (/3/3)oi,-decay of 82 Se and other isotopes as well. 
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calculated in [TTU]. the NEM03 and CUORICINO upper limits read, respectively: |( m }| < (0.50 - 
0.96) eV and |( m )| < (0.25 — 0.43) eV. A very impressive lower limit on the half-life of 136 Xe was 
obtained recently in the EXO-200 experiment [TDU]: T°J 2 ( 136 Xe) > 1.6 x 10 25 yr (90% C.L.). 

The best lower limit on the half-life of 76 Ge, T°J 2 > 1.9 x 10 25 yr (90% C.L.), was found in 
the Heidelberg-Moscow 76 Ge experiment |lllj . It corresponds to the upper limit [1 10] |( m )| < 
(0.20-0.35) eV. A positive (/3/3)o„-decay signal at > 3a, corresponding to Ty 2 = (0.69-4.18) x 10 25 yr 
(99.73% C.L.) and implying |( m )| = (0.1 — 0.9) eV, is claimed to have been observed in j!12j . while a 
later analysis reports evidence for (f3 f3) o u -dec&y at 6<r corresponding to |( m )| = 0.32 ± 0.03 eV 1 13 . 

Most importantly, a large number of projects aim at a sensitivity to |( m )| ~ (0.01 — 0.05) eV |90j : 
CUORE ( 130 Te), GERDA ( 76 Ge), SuperNEMO, EXO ( 136 Xe), MAJORANA ( 76 Ge), MOON ( 100 Mo), 
COBRA ( 116 Cd), XMASS ( 136 Xe), CANDLES ( 48 Ca), KamLAND-Zen ( 136 Xe), SNO+ ( 150 Nd), etc. 
These experiments, in particular, will test the positive result claimed in |113j . 

The existence of significant lower bounds on | ( m ) | in the cases of IH and QD spectra [26] , which 
lie either partially (IH spectrum) or completely (QD spectrum) within the range of sensitivity of the 
next generation of (/3/3)o^-decay experiments, is one of the most important features of the predictions 
of |( m )|. These minimal values are given, up to small corrections, by | Am 2 2 | cos 2#i2 and mo cos 28\2. 
According to the combined analysis of the solar and reactor neutrino data [21J, i) the possibility of 
cos 2#i2 = is excluded at ~ 6a, ii) the best fit value of cos2#i2 is cos2#i2 = 0.39, and iii) at 99.73% 
C.L. one has cos 26\2 ^ 0.28. The quoted results on cos 29\2 together with the range of possible values 
of |Am 2 2 | and mo lead to the conclusion about the existence of significant and robust lower bounds 
on |( m )| in the cases of IH and QD spectrum. At the same time one can always have |( m )| <C 10 -3 
eV in the case of spectrum with normal ordering |49j . As Fig. 1 indicates, |( m )| cannot exceed ~ 5 
meV for NH neutrino mass spectrum. This implies that max{\{ m )|) in the case of NH spectrum is 
considerably smaller than min(\{ m )|) for the IH and QD spectrum. This opens the possibility of 
obtaining information about the type of z/-mass spectrum from a measurement of |( m )\ ^ [261 . 
In particular, a positive result in the future (/3/3)o;,-decay experiments with |( m )| > 0.01 eV would 
imply that the NH spectrum is strongly disfavored (if not excluded). For Am 2 -^^ > 0, such a result 
would mean that the neutrino mass spectrum is with normal ordering, but is not hierarchical. If 
^ m 3i(32) ^ 0, the neutrino mass spectrum would be either IH or QD. Prospective experimental errors 
in the values of oscillation parameters in |( m )| and the sum of neutrino masses, and the uncertainty 



in the relevant NME, can weaken but do not invalidate these results 22 [1TB] [5"T| 154] . 

As Fig. 1 indicates, a measurement of |( m )| k; 0.01 eV would either |49j i) determine a relatively 
narrow interval of possible values of the lightest neutrino mass min(m,j), or ii) would establish an 
upper limit on mm(m,j). If an upper limit on |( m )| is experimentally obtained below 0.01 eV, this 
would lead to a significant upper limit on min(mj). 

The possibility of establishing CP- violation in the lepton sector due to Majorana CPV phases has 
been studied in [49} I52| and in much greater detail in |50t I51j. It was found that it is very challenging: 
it requires quite accurate measurements of |( m )\ (and of mo for QD spectrum), and holds only for a 
limited range of values of the relevant parameters. More specifically |50t I51j. establishing at 2a CP- 
violation associated with Majorana neutrinos in the case of QD spectrum requires for sin 2 00=0.31, 
in particular, a relative experimental error on the measured value of |( m )| and mo smaller than 
15%, a "theoretical uncertainty" F ^ 1.5 in the value of |( m )| due to an imprecise knowledge of the 



22 Encouraging results, in what regards the problem of calculation of the NME, were reported at the MEDEX'll 
Workshop on Matrix Elements for the Double-beta-decay Experiments [114] . For the bounds on |( m )| obtained using 
the current results on the NME see, e.g., [54] , 
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corresponding NME, and value of the relevant Majorana CPV phase 021 typically within the ranges 
of ~ (tt/4 - 3tt/4) and ~ (5tt/4 - 7vr/4). 

The knowledge of NME with sufficiently small uncertainty W\ is crucial for obtaining quantitative 
information on the i/-mixing parameters from a measurement of (/3/3)oi/-decay half-life. The observation 
of a (/3/3)o ! y-decay of one nucleus is likely to lead to the searches and eventually to observation of the 
decay of other nuclei. One can expect that such a progress, in particular, will help to solve completely 
the problem of the sufficiently precise calculation of the nuclear matrix elements for the (/3/3)oi/-decay 



If the future (/3/3)ojy-decay experiments show that |( m )| < 0.01 eV, both the IH and the QD 
spectrum will be ruled out for massive Majorana neutrinos. If in addition it is established in neutrino 
oscillation experiments that the neutrino mass spectrum is with inverted ordering, i.e. that Am^jj < 
0, one would be led to conclude that either the massive neutrinos Uj are Dirac fermions, or that Vj 
are Majorana particles but there are additional contributions to the (/3/3)o^-decay amplitude which 
interfere destructively with that due to the exchange of light massive Majorana neutrinos. The case of 
more than one mechanism generating the (/3/3)oi/-decay was discussed recently in, e.g., |1 1118] . where 
the possibility to identify the mechanisms inducing the decay was also analised. If, however, Am^^^ 
is determined to be positive in neutrino oscillation experiments, the upper limit |( m )\ < 0.01 eV 
would be perfectly compatible with massive Majorana neutrinos possessing NH mass spectrum, or 
mass spectrum with normal ordering but partial hierarchy, and the quest for |( m )| would still be 
open. 

If indeed in the next generation of (/3/3)oi/-decay experiments it is found that |( m )\ < 0.01 eV, 
while the neutrino oscillation experiments show that Am^jj > 0, the next frontier in the searches for 
((3f3)ou— decay would most probably correspond to values of |( m )| ~ 0.001 eV. Taking |( m )| = 0.001 
eV as a reference value, the conditions under which |( m )| in the case of neutrino mass spectrum with 
normal ordering would be guaranteed to satisfy |( m )| ^ 0.001 eV, were investigate in j!02j . In the 
analysis performed in |102j . the specific case of normal hierarchical neutrino mass spectrum, and the 
general case of spectrum with normal ordering, partial hierarchy and values of #13, including the value 



measured in the Daya Bay, RENO, Double Chooz and T2K experiments, eq. (1.4), were considered. 
The ranges of the lightest neutrino mass mi and/or of sin 2 #13, for which |( m )| ^ 0.001 eV were 
derived as well, and the phenomenological implications of such scenarios were discussed. 



1.4 Outlook 

The last 14 years or so witnessed a spectacular experimental progress in the studies of the properties 
of neutrinos. In this period the existence of neutrino oscillations, caused by nonzero neutrino masses 
and neutrino mixing, was established and the parameters which drive the oscillations, were determined 
with a relatively high precision. In spite of these remarkable achievements one has to admit that we 
are still completely ignorant about some of the fundamental aspects of neutrino mixing: the nature - 
Dirac or Majorana, of massive neutrinos, the type of spectrum the neutrino masses obey, the absolute 
scale of neutrino masses, the status of CP symmetry in the lepton sector. Finding out these aspects 
and understanding the origins of the neutrino masses and mixing and the patterns they and possibly 
leptonic CP violation exhibit, requires an extensive and challenging program of research. The main 
goals of such a research program include: 

23 A possible test of the NME calculations is suggested in [49] and is discussed in greater detail in [1 15J . 
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• Determining the nature - Dirac or Majorana, of massive neutrinos Vj. This is of fundamental 
importance for making progress in our understanding of the origin of neutrino masses and mixing 
and of the symmetries governing the lepton sector of particle interactions. 

• Determination of the sign of Arn^jj (Am^) and of the type of neutrino mass spectrum. 

• Determining or obtaining significant constraints on the absolute neutrino mass scale. 

• Determining the status of CP symmetry in the lepton sector. 

• Understanding at a fundamental level the mechanism giving rise to neutrino masses and mixing 
and to L\— non-conservation. This includes understanding the origin of the patterns of neutrino 
mixing and neutrino masses, suggested by the data. Are the observed patterns of z/-mixing and 
of Am|i 31 related to the existence of a new fundamental symmetry of particle interactions? Is 
there any relation between quark mixing and neutrino (lepton) mixing? What is the physical 
origin of CP violation phases in the neutrino mixing matrix Ul Is there any relation (correlation) 
between the (values of) CP violation phases and mixing angles in U? Progress in the theory of 
neutrino mixing might also lead to a better understanding of the mechanism of generation of 
baryon asymmetry of the Universe. 

The successful realization of this research program would be a formidable task and would require 
many years. It already began with the high precision measurement of #13 in the Daya Bay and RENO 



experiments, which showed that sin 2 2^13 has a relatively large value, eq. (1.4). The Double Chooz 
and T2K experiments also found values of sin 2 2#i3, which are different from zero respectively at 3.1a 
and 3.2o" and are compatible with those obtained in the Daya Bay and RENO experiments. These 
results on #13 have far reaching implications. As we have already mentioned or discussed, the measured 
relatively large value of #13 opens up the possibilities, in particular, 

i) for searching for CP violation effects in neutrino oscillation experiments with high intensity 



accelerator neutrino beams, like T2K and NOz^A [H9J 



ii) for determining the sign of Am 2 2 , and thus the type of neutrino mass spectrum, in neutrino 
oscillation experiments with sufficiently long baselines (see, e.g., |23[ 25J). 

A value of sin #13 ^ 0.09 is a necessary condition for a successful "flavoured" leptogenesis with 
hierarchical heavy Majorana neutrinos when the CP violation required for the generation of the matter- 
antimatter asymmetry of the Universe is provided entirely by the Dirac CP violating phase in the 
neutrino mixing matrix |41j. 

With the measurement of #13, the first steps on the long "road" leading to a comprehensive 
understanding of the patterns of neutrino masses and mixing, of their origin and implications, were 
made. The future of neutrino physics is bright. 
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Figure 1.1: The effective Majorana mass |( m }| (including a 2a uncertainty), as a function of min(mj) 
for sin 2 #13 = 0.0236 ±0.0042 [TH] and 5 = 0. The figure is obtained using also the best fit values and 
la errors of Am|i, sin 2 #12, and lAm 2 .-^^! given in Table 1.7 in ref. [BJ. The phases a2i,3i ar e varied 
in the interval [0,7r]. The predictions for the NH, IH and QD spectra are indicated. The red regions 
correspond to at least one of the phases «2i,3i and («3i — o?2i) having a CP violating value, while the 
blue and green areas correspond to 021,31 possessing CP conserving values. (From ref. [B].) 
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